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Abstract 

> 

CN 

' The dynamical properties of an electron-positron-photon plasma created by the 

o ' 

ON ' vacuum polarization process occurring around a charged gravitationally collapsing 

O ' 

core of an initially neutral star are examined within the framework of General Rel- 



ativity and Quantum Field Theory. The Reissner-Nordstrom geometry is assumed 



to apply between the collapsing core and the oppositely charged remnant of the 
star. The appearance of a separatrix at radius R, well outside the asymptotic ap- 
proach to the horizon, is evidenced. The neutral electron-positron-photon plasma 
created at radii r > R self-propels outwards to infinity, following the classical PEM- 
pulse analysis [1,2]. The plasma created at r < R remains trapped and follows the 
gravitational collapse of the core only contributing to the reduction of the electro- 
magnetic energy of the black hole and to the increase of its irreducible mass. This 
phenomenon has consequences for the observational properties of gamma-ray bursts 
and is especially relevant for the theoretical prediction of the temporal and spectral 
structure of the short bursts. 
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The formulation of the physics of the dy ado sphere of an electromagnetic black 
hole (EMBH) has been until now approached by assuming the vacuum polar- 
ization process a la Sauter-Heisenberg-Euler-Schwinger [3,4,5] in the field of 
an already formed Kerr-Newmann [6] or Reissner-Nordstrom black hole [7,8]. 
This acausal approach is certainly valid in order to describe the overall ener- 
getics and the time development of the gamma-ray bursts (GRBs) reaching 
a remarkable agreement between the observations and the theoretical predic- 
tion, in particular with respect to: a) the existence of a proper gamma-ray 
burst (P-GRB) [9], b) the afterglow detailed luminosity function and spec- 
tral properties [10,11,12] and c) the relative intensity of the P-GRB to the 
afterglow [13,10,11]. 

This acausal approach has to be improved by taking into account the causal 
dynamical process of the formation of the dyadosphere as soon as the detailed 
description on timescales of 10~ 4 — 10~ 3 s of the P-GRB are considered. Such 
a description leads to theoretical predictions on the time variability of the 
P-GRB spectra which may become soon testable by a new class of specially 
conceived space missions. 

It is the aim of this letter to report progress in this theoretically challeng- 
ing process which is marked by distinctive and precise quantum and general 
relativistic effects. These new results have been made possible by the recent 
progress in Refs. [14], [8] and especially [15]. There it was demonstrated the 
intrinsic stability of the gravitational amplification of the electromagnetic field 
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at the surface of a charged star core collapsing to an EMBH. The e + e~ plasma 
generated by the vacuum polarization process around the core is entangled in 
the electromagnetic field [16]. The e + e~ pairs do thermalize in an electron- 
positron-photon plasma on a time scale 10 2 — 10 4 times larger than h/m e c 
[15], where c is the speed of light and m e the electron mass. As soon as the 
thermalization has occurred, a dynamical phase of this electrically neutral 
plasma starts following the considerations already discussed in [1,2]. While 
the temporal evolution of the e + e~7 plasma takes place, the gravitationally 
collapsing core moves inwards, giving rise to a further amplified supercriti- 
cal field, which in turn generates a larger amount of e + e~ pairs leading to 
a yet higher temperature in the newly formed e + e~7 plasma. We report, in 
the following, progress in the understanding of this crucial dynamical pro- 
cess: the main difference from the previous treatments is the fact that we do 
not consider an already formed EMBH but we follow the dynamical phase of 
the formation of dyadosphere and of the asymptotic approach to the horizon 
by examining the time varying process at the surface of the gravitationally 
collapsing core. 

The space-time external to the surface of the spherically symmetric collapsing 
core is described by the Reissner-Nordstrom geometry [17] with line element 



with dfl 2 = d6 2 + sin 2 6d(p 2 , a 2 = a 2 {r) = l- 2M/r + Q 2 /r 2 , where M and Q 
are the total energy and charge of the core as measured at infinity. On the core 
surface, which at the time t has radial coordinate r , the electromagnetic field 
strength is 8 = £ (r ) = Qjr\. The equation of core's collapse is (see [14]): 



ds 2 = —a 2 dt 2 + a 



- 2 dr 2 + r 2 dQ 2 , 



(1) 




(2) 
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where H (r ) = — ^f^f and M is the core rest mass. Analytic expressions 
for the solution of Eq. (2) were given in [14] . We here recall that the dyado- 



sphere radius is defined by £ (r ds ) = £ c = m 2 c z /eh [7] as r ds = ^JeQh/mlc 3 , 
where e is the electron charge. In the following we assume that the dyadosphere 
starts to be formed at the instant ids = to (?"ds) = 0. 



Having formulated the core collapse in General Relativity in Eq. (2), in order 
to describe the quantum phenomena, we consider, at each value of r and t , 
a slab of constant coordinate thickness Ar small in comparison with r^s and 
larger than h/m e c 2 . All the results will be shown to be independent on the 
choice of the value of Ar. In each slab the process of vacuum polarization lead- 
ing to e + e~ pair creation is considered. As shown in [15,16] the pairs created 
oscillate [18,19,20,21] with ultrarelativistic velocities and partially annihilate 
into photons; the electric field oscillates around zero and the amplitude of 
such oscillations decreases with a characteristic time of the order of 10 2 — 10 4 
h/m e c 2 . The electric field is effectively screened to the critical value 8 C and the 
pairs thermalize to an e + e~7 plasma. While the average of the electric field £ 
over one oscillation is 0, the average of £ 2 is of the order of £ 2 , therefore the 
energy density in the pairs and photons, as a function of r , is given by [8] 



e (r ) = ± I? 2 (To) - £, 



£ 2 

°c 

Stt 



(3) 



For the number densities of e + e~ pairs and photons at thermal equilibrium 
we have n e + e - ~ n 7 ; correspondingly the equilibrium temperature T , which 
is clearly a function of r and is different for each slab, is such that 



e (T ) = e 7 (T ) + e e + (T ) + e e - (T ) = e , 



(4) 
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with e and n given by Fermi (Bose) integrals (with zero chemical potential): 

e e+e - (To) = ^3 £ ggg^ ^ 2 ^, e 7 (T ) = 4, {kT,f , (5) 
n e+e - (To) = ^ jf n 7 (T ) = ^ (A;T ) 3 , (6) 



where k is the Boltzmann constantFrom the conditions set by Eqs. (4), (5), (6), 
we can now turn to the dynamical evolution of the e + e~7 plasma in each slab. 
We use the covariant conservation of energy momentum and the rate equation 
for the number of pairs in the Reissner-Nordstrom geometry external to the 
star core: 

V a T ab = 0, V a {n e+e -u a ) = av[n 2 e+e -{T)-nl +e J\, (7) 

where T ab = (e + p) u a u b + pg ab is the energy-momentum tensor of the plasma 
with proper energy density e and proper pressure p, u a is the fluid 4— velocity, 
n e + e - is the number of pairs, n e + e - (T) is the equilibrium number of pairs and 
av is the mean of the product of the e + e~ annihilation cross-section and the 
thermal velocity of pairs. We follow closely the treatment which we developed 
for the consideration of a plasma generated in the dyadosphere of an already 
formed EMBH [1,2]. It was shown in [1,2] that the plasma expands as a pair- 
electromagnetic pulse (PEM pulse) of constant thickness in the laboratory 
frame. Since the expansion, hydrodynamical timescale is much larger than the 
pair creation (h/m e c 2 ) and the thermalization (10 2 — 10 4 /i/m e c 2 ) time-scales, 
in each slab the plasma remains at thermal equilibrium in the initial phase of 
the expansion and the right hand side of the rate Eq. (7) is effectively 0, see 
Fig. 24 (second panel) of [11] for details. 

If we denote by £ a the static Killing vector field normalized at unity at spacial 
infinity and by {S t } t the family of space-like hypersurfaces orthogonal to £ a (t 
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being the Killing time) in the Reissner-Nordstrom geometry, from Eqs. (7), the 
following integral conservation laws can be derived (see for instance [22,23]) 



/ i a T ab dL b = E, ( n e+e -u b di: b = N e 



where dS b = a~ 2 £ b r 2 sin OdrdOdcj) is the vector surface element, E the total 
energy and N e + e - the total number of pairs which remain constant in each 
slab. We then have 

\e + p) 1 2 -p]r 2 = (B, n.+^a^r 2 = % s+e -, (9) 

where l£ and 9T e + e - are constants and 



i-«- 4 (i) : 



-1/2 



(10) 



is the Lorentz 7 factor of the slab as measured by static observers. We can 
rewrite Eqs. (8) for each slab as 



(l) 2 =«%, (11) 
( r\ 2 - (t+P) ( n e+e-o \ 2 ( a) 2 P (r\ A ( 12 ) 

\roJ \ £0 / V n e+e- / \ a 0/ e V r / ' V 7 

where pedex refers to quantities evaluated at selected initial times to > 0, 
having assumed r (t ) — r , dr/dt\ t=to = 0, T (t ) = T . 

Eq. (11) is only meaningful when f ro (r) > 0. From the structural analysis of 
such equation it is clearly identifiable a critical radius R such that: 

• for any slab initially located at r > R we have f ro (r) > for any value 
of r > r and f ro (r) < for r < r ; therefore a slab initially located at a 
radial coordinate r > R moves outwards, 

• for any slab initially located at r < R we have f ro (r) > for any value of 
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r+ < r < r and f ro (r) < for r > r ; therefore a slab initially located at a 
radial coordinate r < R moves inwards and is trapped by the gravitational 
field of the collapsing core. 



We define the surface r = R, the dyadosphere trapping surface (DTS). The ra- 
dius R of DTS is generally evaluated by the condition ^ _ = 0. R is so close 

ar r=R 

to the horizon value r + that the initial temperature To satisfies AT~o ^> m e c 2 
and we can obtain for R an analytical expression. Namely the ultrarelativistic 
approximation of all Fermi integrals, Eqs. (5) and (6), is justified and we have 
n e + e - (T) oc T 3 and therefore f ro ~ 1 — (T/T ) 6 (a /a) 2 (r/r ) 4 (r < R). The 
defining equation of R, together with (13), then gives 



R = 2M 



1 + (l - 3Q 2 /AM 2 ) 



1/2' 



> r+. (14) 



In the case of an EMBH with M = 2OM , Q = 0.1M, we compute: 

• the fraction of energy trapped in DTS: 

E=f _ae d£~0.53/ ae rfS; (15) 

Jrj r <r<R Jr + <r<r dB 

• the world-lines of slabs of plasma for selected r in the interval (^R, r ds ) (see 
Fig. 1); 

• the world-lines of slabs of plasma for selected r in the interval (r + , Rj (see 
Fig. 2). 

At time i =t (^R^j when the DTS is formed, the plasma extends over a region 
of space which is almost one order of magnitude larger than the dyadosphere 
and which we define as the effective dyadosphere. The values of the Lorentz 7 
factor, the temperature and e + e" number density in the effective dyadosphere 
are given in Fig. 3. 
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In conclusion we see how the causal description of the dyadosphere formation 
can carry important messages on the time variability and spectral distribution 
of the P-GRB due to quantum effects as well as precise signature of General 
Relativity. 
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Fig. 1. World line of the collapsing charged core (dashed line) as derived from 
Eq. (2) for an EMBH with M = 2OM , Q = 0.1M; world lines of slabs of plasma 
for selected radii ro in the interval r^ s ) . At time i the expanding plasma extends 
over a region which is almost one order of magnitude larger than the dyadosphere. 
The small rectangle in the right bottom is enlarged in Fig. 2. 
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Fig. 2. Enlargement of the small rectangle in the right bottom of Fig. 1. World-lines 
of slabs of plasma for selected radii ro in the interval (r +t R) . 



11 






Fig. 3. Physical parameters in the effective dyadosphere: Lorentz 7 factor, proper 
temperature and proper e + e~ number density as functions at time i for an EMBH 
with M = 2OM and Q = 0.1M. 
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